Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



Supernovae in Deep Hubble Space Telescope 

Galaxy Cluster Fields: Cluster Rates and Field Counts^ 

Avishay Gal-Yam^'^, Dan Maoz^''^, and Keren Sharon^ 

^ Based on observations made with the Hubble Space Telescope, which is operated by AURA, Inc., under NASA contract NAS5-26555. 
^ School of Physics & Astronomy and Wise Observatory, Tel Aviv University, Tel Aviv 69978, Israel; avishay&wise. tau.ac.il 
3 Colton Fellow. 

* Department of Astronomy, Columbia University, 550 W. 120th St., New York, NY 10027, USA; dani&wise. tau.ac.il 
Accepted - . Received - ; 

ABSTRACT 

We have searched for high-redshift supernova (SN) candidates in multiple deep 
Hubble Space Telescope (HST) archival images of nine galaxy-cluster fields. We detect 
six apparent SNe, with 2f.6 < /814 < 28.4 mag. There is roughly 1 SN per deep 
{/814 > 26 mag), doubly- imaged, WFPC2 cluster field. Two SNe are associated with 
cluster galaxies (at redshifts z — 0.18 and z — 0.83), three are probably in galaxies 
not in the clusters (at z = 0.49, z — 0.60, and z = 0.98), and one is at unknown 
z. After accounting for observational efhciencies and uncertainties (statistical and 
systematic) we derive the rate of type-la SNe within the projected central 500/i^q^ 
kpc of rich clusters: R = 0.20t^;^g/iio SNu in 0.18 < z < 0.37 clusters, and R = 
0.4lj;j!|'3gft,gQ SNu in clusters at 0.83 < z < 1.27 (95 per cent confidence interval; 1 
SNu = 1 SN century"-'^ per 10^°Lbo)- Combining the two redshift bins, the mean 
rate is i?z=o.4i = 0.30lg;28/i5o SNu. The upper bounds argue against SNe-Ia being the 
dominant source of the large iron mass observed in the intra-cluster medium. We also 
compare our observed counts of field SNe (i.e., non-cluster SNe of all types) to recent 
model predictions. The observed field count is iV < 1 SN with /814 < 26 mag, and 
1 < < 3 SNe with /814 < 27 mag. These counts are about two times lower than 
some of the predictions. Since the counts at these magnitudes are likely dominated by 
type-ll SNe, our observations may suggest obscuration of distant SNe-II, or a SN-II 
luminosity distribution devoid of a large high-luminosity tail. 

Key words: galaxies: clusters: general - supernovae: general. 



1 INTRODUCTION 

The search for distant supernovae (SNe) has been revolu- 
tionised in the last few years by two cosmology-oriented 
groups, the Supernova Cosmology Project (SCP; Perlmutter 
et al. 1997), and the High-z SN Search Team (Schmidt et 
al. 1998). Both groups have discovered hundreds of distant 
SNe with redshifts up to 2 = 1, and beyond. Analysis of 
these data shows evidence for an accelerating Universe, pos- 
sibly driven by a positive cosmological constant (Riess et al. 
1998; Perlmutter et al. 1999). Two likely high-z SNe have 
also been discovered in a repeated observation with the Hub- 
ble Space Telescope {HST) of the Hubble Deep Field (HDF; 
Gilliland, Nugent, & Phillips 1999). One of these SNe has 
been shown by Riess et al. (2001) to be a type-la a,t z — 1.7, 
at an epoch before the transition from deceleration to ac- 
celeration. A third possible SN candidate was detected in 
the original HDF data through the analysis of the time de- 



pendent flux variation of objects in this field (Mannucci & 
Ferrara 1999). These observations demonstrate that deep 
HST observations are an efficient way to find high— z SNe, 
with the high angular resolution and low sky background 
compensating for the small field of view. 

Apart from the use of type-la SNe as standard can- 
dles, it is important to find and characterise distant SNe of 
all types. SN rates as a function of redshift probe the star 
formation history of the Universe, the physical mechanisms 
leading to SN-Ia explosions, and the cosmological parame- 
ters (Jorgensen et al. 1997; Sadat et al. 1998; Ruiz-Lapuente 
& Canal 1998; Madau, Delia Valle, & Panagia 1998; Yun- 
gelson & Livio 2000). Furthermore, even in the absence of 
redshift information, SN number counts as a function of lim- 
iting magnitude can provide valuable constraints (Dahlen & 
Fransson 1999, hereafter DF; Sullivan et al. 2000, hereafter 
S2000). 
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There are incentives to search for high— 2; SNe not only 
in "blank" fields, but also in the fields of rich galaxy clusters. 
The lensing magnification of background sources by the clus- 
ter potentials enables detection of distant and intrinsically 
dim sources, that would be too faint to study otherwise. In- 
deed, the natural telescopes provided by rich clusters have 
been used successfully for deep studies of high— z galaxies 
in the X-ray (Crawford et al. 2001), optical (e.g., Frye & 
Broadhurst 1998, Pettini et al. 2000), IR (e.g., Lemonon et 
al. 1998), and submm (Small et al. 1998) regimes. Follow- 
ing the same reasoning, S2000 have considered the prospects 
of searching for high— 2 SNe lensed by galaxy clusters. The 
most distant SNe could be found, with implications both for 
the cosmological applications and for our understanding of 
the nature of SNe. 

Lensed high— z SNe behind galaxy clusters could also 
probe the properties of the cluster lenses. For instance, 
a magnified type-la SN would provide a direct measure- 
ment of the cluster magnification, lifting the so called 
"mass-sheet degeneracy" encountered when deducing clus- 
ter masses from weak-lensing shear maps (Kolatt & Bartel- 
mann 1998). A SN in a host galaxy that is lensed by a cluster 
into a giant arc may be multiply imaged into three images 
with time delays of weeks to months (Kovner & Paczynski 
1988). SNe in galaxies that are nearly lensed into optical 
Einstein rings, such as the CL0024-(-16 arc system (Tyson, 
Kochanski, & dell' Antonio 1998), may also have short time 
lags (months to years) between images of the SNe in distinct 
galaxy images, providing constraints on the cluster mass dis- 
tribution and on cosmological parameters. 

Finally, knowledge of the rate of SNe in the galaxy clus- 
ters, as a function of SN type and cluster redshift, could illu- 
minate several puzzles concerning the intra-cluster medium 
(ICM). One such puzzle is the total observed mass of iron, 
which is at least several times larger than that expected 
from SNe, based on the present-day stellar masses, and as- 
suming a standard stellar initial mass function (IMF; e.g., 
Brighenti and Mathews 1998; Loewenstein 2000). Another 
is the energy budget of the ICM gas and the "entropy fioor" 
observed in clusters, which suggest a non-gravitational en- 
ergy source to the ICM, again several times larger than the 
expected energy input from SNe (e.g., Lloyd-Davies, Pon- 
man, & Cannon 2000; Tozzi & Norman 2001; Brighenti & 
Mathews 2001). Proposed solutions to these problems have 
included an IMF skewed toward high-mass stars (so that 
a large number of iron-enriching core-collapse SNe are pro- 
duced per present-day unit stellar luminosity), or a dom- 
inant role for SNe-Ia in the ICM enrichment. The latter 
option would require early-type galaxies in clusters to have 
SN-Ia rates much higher than observed in nearby field el- 
lipticals. Constraints on these scenarios can be obtained by 
directly measuring the SN rate in clusters at various red- 
shifts. 

Observationally, SN searches in galaxy clusters were pi- 
oneered in the late 1980's by Norgaard-Nielsen et al. (1989), 
resulting in the first detection of a z = 0.31 SN in the galaxy 
cluster AC118. More recently, low-redshift cluster samples 
have been monitored for SNe by the Mount Stromlo Abell 
Cluster SN Search (Reiss et al. 1998) and by the Wise Obser- 
vatory Optical Transients Search (Gal- Yam & Maoz 1999, 
2001). To extend these efforts to higher— z clusters, we have 
carried out a search for SNe in duplicate deep cluster images 



in the HST archive. In this paper, we describe our discov- 
ery of six SNe in this data set. We then use these SNe to 
derive the SN-Ia rate in clusters and to test the theoretical 
predictions for the non-cluster SN number counts in such 
fields. Where applicable, we assume a flat cosmology with 
parameters i},n = 0.3 and Ha = 0.7, and use hso to denote 
the Hubble parameter in units of 50 km s~^ Mpc~^. 



2 ARCHIVAL TARGET SELECTION AND 
ANALYSIS 

We have searched the HST archive for duplicate WFPC2 
observations of galaxy cluster fields. Clusters were identified 
as such by the target name, or, when the names were not re- 
vealing, by searching the NED database for known clusters 
near the field coordinates. The details of our final cluster 
sample are given in Table 1. For two WFPC2 observations 
to be considered duplicates, we required that they had been 
obtained through the same filter, and that the angular dis- 
tance between their field centres was less than 75 arcsec. The 
requirement for identical filters is due to the need for accu- 
rate image subtraction when searching for faint transients. 
The angular separation criterion is driven by the peculiar 
shape of the WFPC2 field, coupled with the random HST 
orientations, which can result in observations that are nom- 
inally duplicate, but have little overlap area. 

To optimise the search for SNe, especially at high z, 
where cosmological time-dilation is important, we required 
a minimal time separation between data sets of 30 days. The 
implications of this criterion are accounted for in the analy- 
sis of our results. Finally, we have limited our search to deep 
observations, with total exposure times of at least 2000 s in 
one of the broad-band HST filters. In order to facilitate rea- 
sonable rejection of cosmic rays and hot pixels, we required 
that at least one of the duplicate data sets consist of at 
least three sub-exposures. When both duplicate data sets 
had three or more sub-exposures, we searched for candidate 
variable objects in both. In cases where only one of the data 
sets was suitably split, the other set was used only as a ref- 
erence, since the possibility of cosmic ray events mimicking 
transient objects could not be ruled out. We have found nine 
galaxy cluster fields with observations in the HST archive 
satisfying these criteria. The overlap between images from 
different epochs is partial, so the total effective area useful 
for a bi-directional search is 30.8 arcmin^, or 6.6 times the 
total area of WFPC2, excluding the PCI CCD (see Table 
!)• 

Sets of WFPC2 images at a given epoch were com- 
bined using software based mainly on the DITHER pack- 
age (Fruchter & Hook, 1997) within IRAF 0. A flexible, 
high-fldelity algorithm was used to reject cosmic-rays and 
hot/cold pixels, even in stacks with as few as three sub- 
exposures, allowing their use in this project. The PCl-CCD 
data were not used in the analysis, due to the lower sen- 
sitivity, the relatively small addition of area, and the need 

* IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, wliicli are op- 
erated by AURA, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Galaxy Cluster Fields with Duplicate i75T Observations 



Cluster 


z 


Epoch 


Exposure 


Sub- 


Filter 


Program 


Overlap Area" 


Number of 






(UT Start) 


[ks] 


frames 




ID 


[arcmin^] 


Apparent SNe 


Abell 2218 


0.18 


Mar 11 1999 


8.4 


12 


F606W 


7349 


4.69 


1 






Jan 13 2000 


10.0 


10 


F606W 


8500 


4.69 


- 


Abell 1689 


0.18 


Jun 27 1996 


2.0 


2 


F814W 


6004 





- 






Jun 19 1997 


4.0 


4 


F814W 


6004 


4.69 


- 






Jul 07 1997 


4.6 


4 


F814W 


5993 


1.56 


- 






Feb 24 2000 


4.0 


4 


F814W 


8546 


1.56 


- 


AC114 


0.31 


Jan 06 1996 


16.6 


6 


F702W 


5935 


1.28 


- 






Oct 26 1997 


15.6 


6 


F702W 


7201 


1.28 


- 


MS1512. 4+3647 


0.37 


Mar 19 1996 


6.3 


9 


F555W 


6003 


2.30 


- 






Jun 20 1997 


10.4 


8 


F555W 


6832 


2.30 




MS2053.7— 0449 


0.58 


Oct 23 1995 


2.6 


2 


F814W 


5991 











Sep 25 1998 


3.2 


3 


F814W 


6745 


4.42 




MS1054.4-0321 


0.83 


Mar 13 1996 


15.6 


6 


F814W 


5987 


4.69 


3 






May 29 1998 


6.6 


6 


F814W 


7372 


4.69 




CL1604+4304 


0.89 


Feb 25 1994 


32.0 


16 


F814W 


5234 


3.97 


1 






Feb 28 1995 


32.0 


16 


F814W 


5234 


3.97 


1 


RXJ0848.6+4453 


1.27 


Mar 01 1999 


13.9 


5 


F814W 


6812 


4.69 








Apr 15 1999 


13.9 


5 


F814W 


6812 


4.69 




CL1645+46 


> 2.8 


May 7 1997 


6.0 


12 


F814W 


6598 


3.42 








Sep 30 1999 


5.3 


4 


F814W 


7342 


3.42 





Notes: 

" The overlap area is set to zero for data sets with fewer then three sub-exposures. These data sets are used 
only as references for searching for SNe at other epochs (§ 3.3.3). 



to resample the data in order to compare to WF-CCD data 
from other epochs. 

The combined images of a cluster from two available 
epochs were registered, using a general geometrical trans- 
formation, calculated from the positions of ~ 20 compact 
galaxies or stars, visible at both epochs. The images were 
scaled according to their exposure times and subtracted, 
forming a difference image. The stable HST point-spread 
function (PSF) allows for good cancellation of almost all 
non-variable objects, as shown in an example in Fig. 1. 
Some low-level residuals remain in the difference images 
near bright stars and at the nuclei of the brightest galax- 
ies. These residuals are unlikely to mimic real objects, since 
their suspect origin is evident from the images. However, the 
residuals near galactic nuclei limit the detection of nuclear 
variability in these galaxies, either SNe close to the nucleus 
or active galactic nuclei (AGN). This is taken into account 
by our detection efficiency simulations, described in § 3.3.2. 

The difference images were methodically scanned by 
eye, and all residual objects checked and classified. After 
the elimination of obvious subtraction artefacts, the loca- 
tion of each residual object was checked in each of the sub- 
exposures from that epoch. For a candidate to be considered 
secure, we required that it have a stellar PSF, and appear in 
at least three independent subsets of the data, each consist- 
ing of the sum of individual exposures in which its location 
is not obviously affected by cosmic rays, with non-variable 
flux levels, to within errors. This criterion eliminates residual 
events consisting of cosmic rays that are coincident in several 
sub-exposures, as well as chip defects and rare, high ampli- 
tude, background fluctuations, since in these cases most of 
the candidate's flux in the final combined image comes from 
only some of the available sub-exposures. 



Care was taken not to reject real events affected by cos- 
mic ray hits, even in data sets with few sub-exposures. Real 
events could, in principle, be lost in three-image data sets if 
cosmic rays happened to hit the object's location on two of 
the three sub-exposures. However, given that only 1-2 per 
cent of the pixels in each image are affected by cosmic rays, 
this has a probability of ~ 10""*, and much lower in data 
sets with more than three sub-exposures. Candidates that 
passed these criteria were considered real transient or vari- 
able sources, most likely high— 2 SNe (see §3.2). In practice, 
all but one of the transient sources we found can be identi- 
fied on each individual exposure, and not only in summed 
subsets. Figure 2 shows an example. 

The fiux limit at each epoch depends on the length of 
the exposures and the filter used, and therefore varies be- 
tween our sample clusters, but in all sets the limit for > 50 
per cent detection is better than /814 — 26 maaU 



3 RESULTS 

3.1 Apparent supernovae 

We have detected six apparent SNe in the cluster fields we 
have searched, including the rediscovery of SN 1996cl, which 
was previously known. For each SN, Table 2 lists the details, 
and Fig. 3 shows a section of the images at two epochs. 
Among the six SNe, five are projected within 2.5 half-light 

t We report our observations using the AB magnitude system, 
in which the magnitude in a band b centered at wavelength A 
[nm] is bx = —48.6 — 2.5 X log/^, where fi, is the flux density at 
frequency u in erg cm~^ Hz~^. 
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Figure 1. Example of image subtraction; SN 1996cl, in a galaxy 
in cluster MS1054.4-0321, in a 8'.'5 x 3'.'6 section of the difference 
image (lower panel) created by subtracting the data obtained in 
May 1998 (upper panel) from the data of the same field taken 
during March 1996 (middle panel). Note the clean subtraction of 
the galaxies, resulting from the stable HST point-spread function. 

radii of galaxies that are likely their hosts, while the sixth 
has no obvious host (but see § 3.1.2). We have been able to 
determine the rodshift and typos of the host galaxies of five 
events from the literature (see below). As noted in proof by 
S2000, the object they previously considered a SN candidate, 
near the core of AC114, is not variable. Indeed, Campusano 
et al. (2001) show that this object is likely one of the images 
of a multiply-lensed background galaxy (marked A4 in their 
Fig. 1). We briefly discuss each SN below. 

3.1.1 SN 1 999gv in the field of Abell 2218 

During March 1999, a large mosaic (22 positions) of WFPC2 
images, covering the entire field of Abell 2218, was obtained 
through the F606W filter. The main goal of this program 
(PI Squires) was to create weak Icnsiug shear maps, and 
to compare the resulting mass distribution to that derived 
via strong lensing. X-rays, and the Sunyaev-Zeldovich ef- 
fect. In January 2000, following the Space Shuttle servicing 
mission 3a, the cluster was re-imaged through the F450W, 
F606W, and F814W filters as an "early release" observation, 
intended to demonstrate the telescope performance. Com- 
parison between the 8.4-ks F606W image from March 1999 
and the 10.0-ks F606W image from .January 2000, reveals a 
bright unresolved object in the 1999 exposure (Gal- Yam & 
Maoz 2000; see Fig. 3). The host galaxy is catalogued as No. 
149 in Le Borgne et al. (1992), where it is classified as an 




Figure 2. Example of candidate verification; 2'.'1 X 2'.'1 sections 
of the individual WFPC2 exposures at the location of SN 1995bf 
in CL 1604-1-4304. The apparent SN is visible in all sub-exposures 
where the location is not affected by cosmic-ray hits, with a stellar 
PSF and constant flux, to within the expected errors. 

elliptical galaxy with B = 19.9 mag and redshift z = 0.1753. 
In view of the brightness of this SN, and its association with 
an elliptical cluster member, it is most likely a type-la SN 
in the cluster. Assuming a typical SN-Ia, the apparent mag- 
nitude suggests a SN age of ~ 25 days past maximum light. 

3.1.2 SNe 1996d, 1996cp, and 1996cq in the field of 

MS1054.4- 0321 

The core of the distant galaxy cluster MS1054.4— 0321 {z = 
0.83) was observed with the F814W filter during March 
1996, in a program to study distant. X-ray selected clus- 
ters (PI Donahue). It was re-observed in May 1998, as part 
of a shallower, wide field mosaic (PI Franx), designed for a 
morphological large-scale study. A comparison of the images 
from the two epochs reveals three apparent SNe (Gal- Yam, 
Sharon, & Maoz 2001). 

The brightest of the three is SN 1996cl, originally dis- 
covered by the SCP in ground-based data obtained on March 
18, 1996 (Perlmutter et al. 1996), and spoctroscopically con- 
firmed as a type-la SN at the cluster rodshift. In the HST 
observation, obtained 5 days prior to that of Perlmutter et 
al., the SN is detected with 78i4 = 23.71 mag, the magni- 
tude expected for SN-Ia some 10 days (~ 5.5 days in the SN 
rest frame) before peak magnitude. 

A second transient object in this field, SN 1996cq, is 
detected at /gi4 = 25.61 mag, without any obvious host 
galaxy. This object is superposed on a faint fuzz near the 
image detection limit that may be a low-surface-brightness 
galaxy associated with it. 

The third transient source is SN 1996cp, with Isi4 = 
26.7 mag. This object is superposed on a small knot of emis- 
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Figure 3. Sections of tiie images, at two epochs, for each of the six apparent SNe. The scales shown in the upper-left-hand corners 
correspond to 1". 



Table 2. Apparent Siipcriiovac 



SN 


Cluster 


Magnitude 


a 


5 


Offset 


Host 




Field 




(2000) 


(2000) 


from Host 


2 


1999gv 


Abell 2218 


VfioG = 21.64 ±0.03 


16'"36"02.2= 


+66°12'34" 


0'.'8 S 2'.' 8 W 


0.175 


1996cl 


MS1054.4-0321 


/814 = 23.71 ± 0.04 


10'" 56™ 59.1= 


-3°37'36" 


O'.'l S C.'OS w 


0.827 


1996cq 


MS1054.4-0321 


/814 = 25.61 ± 0.05 


10'"57"02.2= 


-3°37'34" 




? 


1996cp 


MS1054.4-0321 


/si4 = 26.7 ±0.2 


10'"57™03.5'' 


-3°37'41" 


0'.'9 E 


0.596'" 


1995bf 


CL 1604+4304 


/814 = 26.14 ±0.06 




+43°04'44" 


2'.'9 N 0'.'6 W 


0.985 


1994ao 


CL1604+4304 


/814 = 28.4 ± 0.5 


Ig/j04m23.3s 


+43°03'56" 


1'.'4 S I'.'l W 


0.496 



Notes; 

" Other host galaxies possible, see text. 



sion, 0.9 arcsec due east from the nucleus of a blue, late type 
galeixy with z = 0.60 (Milvang- Jensen et al. 2002). However, 
it could also be associated with a larger spiral galaxy that is 
found 4.1 arcsec to the south-east, catalogued as No. 1403 
by van Dokkum ct al. (2000), who classify it as a cluster 
member with z = 0.81. Since the emission knot which hosts 
SN 1996cp is closer to the z = 0.60 galaxy and has similar 
.^814 — Veoe colours, we adopt a redshift of z = 0.60 for this 
SN. However, we cannot at this time exclude the possibil- 



ity that the SN is associated with the nearby cluster galaxy 
(No. 1403) or that the emission knot is a separate galaxy 
altogether. Derivation of photometric redshifts for the emis- 
sion knot hosting SN 1996cp from multi-color images of this 
cluster may shed more light on this issue. 
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3.1.3 SNe 1994ao and 1995bf m the field of CL1604+4304 

The cluster CL1604+4304 (z = 0.89) was the target of a 
deep observation in 1994 (32 ks through the F814W filter; 
PI Westphal). However, a pointing error placed the target 
~ 20 arcsec from the designated position. The cluster was 
therefore re-imaged one year later, again for 32 ks. Subtrac- 
tion of these two deep images reveals two SNe (Gal- Yam 
et al. 2001). The first, SN 1995bf, with Igu = 26.14 mag, 
is found near a putative host, catalogued by Lubin et al. 
(1998) as No. 019, an Sb galaxy with z = 0.985, i.e., behind 
the cluster. 

The limiting magnitude for SN detection in these par- 
ticular images, derived from our detection-efficiency simula- 
tions (see §3.3.2), is Isia = 29 mag, if the requirement that 
the SN appear in individual sets of sub-exposures is relaxed, 
to the effect that candidates need be detected with a stellar 
PSF and constant flux in only two independent subsets of 
the data. This analysis reveals a second, very faint, appar- 
ent SN in this chistcr. SN 1994ao has /gi4 = 28.4 mag, near 
the detection limit, and is the faintest SN detected to date 
(~ 1.5 mag fainter than SN 1997ff; Gilliland et al. 1999; 
Riess et al. 2001). This object is near an SO gala:x;y with 
z = 0.496, in the cluster foreground (No. 028 of Lubin et 
al. 1998). If the event is indeed associated with this host, 
it could imply it is a SN la, as core-collapse SNe (types II, 
lb, Ic) are generally not found in SO galaxies. However, the 
absence of core-collapse SNe in SO galaxies has only been 
established in low redshift (z < 0.1) environments. If this 
is indeed a typo-Ia SN, then the low luminosity implies it 
is many months past maximum. We report this event as an 
example of the capability of HST to find very faint SNe, but 
do not include it in our quantitative analysis of the results, 
below. 



3.2 Supernova alternatives 

Apart from SN 1996cl, all the apparent SNe we discuss lack 
spectroscopic confirmation. We now consider alternative ex- 
planations, other than SNe, for the observed objects, and 
argue that all are unlikely. 

Solar system objects — Small bodies in the Solar 
System (e.g., Kuiper Belt objects - KBOs; Trujillo, Jewitt 
& Luu 2001, and references therein) can potentially mimic 
transients. To pose as one of our apparent SNe, such an 
object must be faint (/8i4 > 22 mag) and slowly moving. 
However, the combination of relatively long observing peri- 
ods (hours to days) with the HST resolution allows us to 
set firm upper limits on the proper motions of the detected 
transient events. For the relevant data set with the shortest 
observing period (MS1054.4— 0321), the upper limit on the 
motion of the transients is /u < 0.12"day~^. For all three 
fields in which we have found SN candidates, the expected 
parallaxes of a KBO at 50 AU, in the field directions and at 
the times of observation, as a result of the Earth's orbital 
motion, are > 70"day~^. This is much larger than our upper 
limit. We note that the proper motion of a KBO at several 
tens of AU from the Sun cannot be large enough to offset 
the effect of the Earth's parallax, so that KBOs, regardless 
of their exact orbital parameters, cannot account for any of 
the events we discuss. 



Variable Galactic stars — The fainter (/8i4 ~ 
26 mag) among the transient events we have discovered are 
generally too faint to be plausible variable Galactic objects. 
Such objects would need to vary by at least ^ 1 mag in 
order to fall below the detection limit on one of the epochs. 
A chromospherically active I = 27 mag M4 dwarf (which, 
say, flared into our detection range), would be at a distance 
of 40 kpc, far out in the low-stellar-density halo, and hence 
highly unlikely to turn up in such a small field. Cataclysmic 
variables in quiescence will always be brighter at minimum 
hght than a single M-dwarf. Our only SN candidate brighter 
than 7814 = 25.5 that lacks spectroscopic confirmation, SN 
1999gv in Abell 2218, coxild be a bright flare from a Galac- 
tic G- or K-dwarf. However, the lack of any evidence for 
such flares at this position in numerous previous images of 
Abell 2218, obtained over several years through various tele- 
scopes (e.g., Gal- Yam &: Maoz 2000), makes this possibility 
unlikely. 

Vciriable stars in other galcLxies — If the transient 
events we have found are associated with the galaxies near 

to which they are projected, then they are generally too lu- 
minous to bo anything but SNe. Of the six events we discuss, 
the redshifts of the probable hosts of five events are known 
(Table 2). The most luminous novae have absolute magni- 
tudes Mi = -11 (Delia VaUe 1991; Warner 1995) and thus 
would be far too faint to bo detected at the redshifts of these 
galaxies, falling below I ~ 28.5 mag at z ~ 0.17. The host 
galaxy of SN 1996cq is only marginally detected in the deep 
HST image. If this event is a nova, the implied luminosity 
of the host becomes implausibly low. Mi > —8.7 mag. 

Active galactic nuclei — None of the events we dis- 
cuss are located in the nuclei of the putative host galaxies, 
and hence are not Seyfert-like AGN. The candidates could, 
in principle, be more distant, highly variable quasars in the 
nuclei of undetected galaxies, projected by chance within 
< 2.5 half-light radii of the galaxies that are seen in five of six 
cases. By performing aperture photometry for all the galax- 
ies in each image, we have determined the galaxies' half-light 
radii. We find that, even in the most crowded fields, < 5 per 
cent of the total area we have searched is within 2.5 half- 
light radii of a galaxy. For SN 1996cq, whose candidate host 
is barely detected, the quasar option cannot be ruled out as 
firmly. 

Gamma-ray burst afterglows — Recent observations 

of gamma-ray burst optical afterglows suggest that these ex- 
plosions are not isotropic (e.g., Frail et al. 2001, and refer- 
ences therein). If the gamma rays are more beamed than the 
optical emission, there may be many optical afterglows that 
have no corresponding gamma-ray burst (Rhoads 1997). 
Such events could, in principle, explain some of the tran- 
sients we have discovered. However, recent work by Dalai, 
Gricst, & Pruet (2001) shows that even a continuous moni- 
toring campaign reaching a limiting magnitude of i? = 27 is 
expected to yield < 0.15 events deg^^ yr^^. Our search has 
visited a total area < 0.01 deg^, only twice. 

We conclude that most or all of the events we have dis- 
covered are likely bona-fide SNe. Nevertheless, in our deriva- 
tion of cluster SN-Ia rates and field counts, the estimated 
uncertainties account for the more ambiguous identifications 
of several of the candidates. In our comparisons of field 
counts to some theoretical predictions (§3.3.3), if any of the 
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events were not SNe, the discrepancy between observations 
and some models would only grow. 

3.3 Field supernova counts 

3.3.1 Observed counts 

Of the six apparent SNe we have detected, three have Isi4 < 
26 mag, five have /si4 < 27 mag, and one has Isi4 ~ 28.4 
mag. SN 1996cl in MS1054.4-0321 is a spectroscopically 
verified cluster event, and SN 1999gv in Abell 2218 is also 
most likely associated with a cluster galaxy. Thus, our upper 
limits on the detected field (i.e., non-cluster) SN counts, is 
one event with 7314 < 26 and three events with Isi4 < 27. 
These observational results can bo compared directly with 
theoretical predictions for the number of SNe expected in 
deep WFPC2 images. 

The data in our sample were obtained with a variety of 
filters and exposure times. Converting our results to absolute 
number counts (i.e., SNe per unit sky area, visible at some 
moment to a given flux limit) is further complicated by a 
dependence on poorly known parameters for core-collapse 
SNe. Instead, we take the published theoretical predictions, 
fold them through our observational parameters, apply the 
model assumptions self-consistently, and derive the expected 
numbers of SNe in our particular sample. We then test the 
models by direct comparison of the expected numbers to the 
observations. 

3.3.2 Detection efficiencies 

We have carried out simulations in order to estimate our SN 
detection efficiencies. Several tens of artificial point sources 
("fake" SNe) with known fiux and the proper amount of 
Poisson noise were blindly added to the data for each field. 
The fake SNe were distributed such that the probability that 
a galaxy hosts a fake SN is proportional to the galaxy bright- 
ness. The separation of a SN from the galaxy centroid was 
chosen from a Gaussian distribution with a equal to 0.85 
times the galaxy's half-light radius. Thus, the fake-SN dis- 
tribution follows the galactic light. Although in reality the 
SNe should follow the galactic luminosity, flux is a reason- 
able surrogate; the cluster galaxy population in a given field 
is all at a similar distance, and most of the field galaxies 
are likely also in a limited redshift range. Our procedure is 
also conservative, in the sense that it places the fake SNe on 
too-high galactic backgrounds, leading to an underestimate 
of the true efficiency, since cosmological surface-brightness 
dimming increases the contrast between a SN and its host- 
galaxy background at high rodshifts. In principle, one might 
place some fake SNe at blank sky positions, to represent 
SNe in host galaxies that are below the detection limit. We 
find that our detection efficiency for such objects is sig- 
nificantly higher than for objects placed on galactic back- 
grounds. However, having no concrete estimate for the frac- 
tion of "hostless" SNe from the total population, and in line 
with our conservative approach, we did not include such 
events in the final simulations. 

The simulated data wore reduced and searched for SNe, 
like the real data, and the fraction of recovered objects 
noted. Examples of the efficiency curves for the shallow- 
est (Abell 1689) and deepest (CL1604-I-4304) data sets are 



Table 3. Sample Detection-Efficiency Values 



Magnitude Efficiency 



^814 


Abell 1689 


CL1604-I-4304 


23 


0.95 


1.0 


24 


0.90 


1.0 


25 


0.80 


1.0 


26 


0.45 


0.85 


27 





0.60 


28 





0.50 


29 









given in Table 3, and the relevant efficiency values for all 
fields are shown in Table 4. 

Six of the nine clusters we discuss were observed 
through the F814W filter, which resembles the Kron- 
Cousins I filter. Model predictions for SN number counts 
in deep HST fields have been calculated for the I band (DF; 
S2000; see §3.3.3), and for these fields the comparison is 
straightforward. However, for three of the cluster fields, the 
duplicate observations were obtained through the F555W 
(Johnson V), F606W ("wide F"), and F702W ("wide i?") 
filters. 

A quantitative comparison of the observed V-band or 
_R-band number counts with the expected J-band counts re- 
quires knowledge of the characteristic colours of the high— a 
SN population. While this information is fairly well known 
for type-la SNe (Perlmutter et al. 1999, Schmidt et al. 1998), 
beyond 25-26 mag, the counts are expected to be dominated 
by core-collapse SNe (e.g., DF; S2000). Broad-band photom- 
etry and spectroscopy of core-collapse SNe with z > 0.5, as 
well as UV spectra of local core-collapse SNe, for deriving 
K-corrections, are scarce (DF; Gilliland et al. 1999; S2000). 
As the colours of SNe depend strongly on the SN redshift, 
type, age, and extinction, the colour terms are highly model 
dependent. The relative fractions of the various sub- types 
of core-collapse SNe (II-P, II-L, Iln, lb, and Ic) are observa- 
tionally unknown at high-z, and must be assumed based on 
sparse data from local observations. With these limitation 
in mind, we pursue two alternative analyses. 

First, following DF and S2000, we assume that the SN 
population probed by the data is dominated by type-II SNe 
aX z ^ 0.7, about half of which are type II-P events. (The 
actual fraction of the various SN sub-types is controversial; 
we merely mate the same assumptions as DF and S2000 
in order to compare their models and our observations self- 
consistently.) Based on UV spectra of local type-II SNe, DF 
assume that type II-P SNe show a strong UV deficit as they 
enter the plateau phase, some 20 days after maximum light, 
while type-II-L and -Iln SNe do not have this effect. At high 
z, the type II-P SNe will therefore have red V — I and R — I 
colours. We estimate that about half of the type II-P SNe are 
observed close to maximum light, when V — I and R — I are 
close to zero (Fig. 10 of DF), and apply no colour correction 
to the efficiency for this fraction of the population. For the 
remaining type II-P SNe, observed at rest-frame ages t >20 
days, we adopt the colours calculated by DF, based on the 
SN models by Eastman et al. (1994). Converting the colours 
from Fig. 10 of DF to the magnitude system we use, we 
obtain Ryo2 — Isu = 0.75 and — Isu > 2. Since the 
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F606W filter extends further into tfie red tlian F555W, and 
lias significantly higher sensitivity, we assume Veoe — ^8i4 ~ 
1.5. The detection efficiency, ri, that applies to the half of the 
type II-P SN population with ages t > 20 days (1/4 of the 
total SN population) is the efficiency found for the observed 
(non-7814) band, at a magnitude obtained by adding the 
colour term. For example, in the Abell 2218 data set, the 
corrected efficiency at /8i4 = 26.4 mag is [0.75 x »7(V606 = 
26.4)] + [0.25 X r;(V606 = 26.4 + 1.5)] = [0.75 x 0.63] + [0.25 x 
0] = 0.47. 

Alternatively, in order to compare our observed SN 
number counts with the predictions in the least model- 
dcpcndcnt way, we consider only the sub-sample of six fields 
observed through the F814W filter, avoiding the problem of 
colour corrections altogether. This approach leads to slightly 
smaller number statistics, but the conclusions, presented be- 
low, remain unchanged. 

In all but one of the data sets we have analyzed, the 
interval between epochs is ten months or longer, ample time 
for even the slowest SNe to vary significantly, even when 
time dilation at z ~ 1 is taken into account. The two epochs 
we compare in the CEise of RXJ0848. 6+4453, however, are 
separated by only 44 days, equivalent to 22 rest-frame days 
for a, z = 1 SN. Slowly- varying SNe Iln and SNe II-P during 
their plateau phase could become hard to detect in this data 
set. To correct for this possible bias, we again assume that 
about 1/2 of the high— « SNe we expect to discover are SNe 
II-P. These are expected to be in the slowly varying phase of 
their light curve for about 1/3 of the time they are observ- 
able. We therefore estimate that 1/6 of the expected high— a 
SNe would not be detected, and apply a correction of 5/6 to 
the core-collapse SN detection efficiency in this data set. One 
also needs to correct for the possibility that SNe observed 
on the rise in the first epoch would have approximately the 
same magnitude (within our observational error) when ob- 
served in the declining phase 44 days later, and thus avoid 
detection. Using simulated variable objects superposed on 
the actual images, we find that 7 = 26 objects can vary by 
at most 0.4 mag and still avoid detection. Applying this to 
typical light-curves for the various SN types, we estimate 
the about 15 per cent of the entire SN population might be 
lost. We therefore apply a further correction to our efficiency, 
with the final correction factor being (5/6) x 0.85 = 0.7 . 

3.3.3 Comparison with theoretical predictions 

S2000 have considered in detail the exact observational route 
we have taken - a search for SNe in deep, duplicate HST 
observations of galaxy cluster fields. Using models for the 
expected SN rates, and taking into account various mod- 
els for the lensing effects of the clusters, S2000 have cal- 
culated the number of SNe expected when comparing two 
images of the same cluster from different epochs. They find 
that the SN number counts depend weakly on the lensing 
and star- formation models they use (Table 1 of S2000). The 
background SN number counts are weakly affected by lens- 
ing because the added depth given by the magnification is 
countered by the reduced effective area that is searched, 
due to that same magnification (Poiciaiii & Madau 1999; 
S2000). We can therefore adopt representative values from 
S2000 for the number of SNe expected per WFPC2 field for 
a two-epoch visit. Converting the Vegarbased magnitudes of 



S2000 (M. Sullivan, private communication) to AB magni- 
tudes, their predicted counts arc Np = 1.5 SNe per WFPC2 
field to 7814 < 26.4 mag, and iVp = 3 SNe to 78i4 < 27.4 
mag, in a bi-directional search (i.e., including SNe discov- 
ered in both epochs). Note that these numbers arc for field 
SNe only, as S2000 exclude cluster events from their calcula- 
tions. S2000 assume 5.71 arcmin^ for the area of a WFPC2 
field, while the correct value for the WF CCDs only is 4.69 
arcmin^, so we scale their predictions by 0.82. 

In similar work, DF have also calculated the num- 
bers and properties of high— 2; SNe. Although they address 
mainly the expected number of SNe detectable with the 
Next Generation Space Telescope, in their Table 2 they also 
give the expected number of SNe in an HST WFPC2 field 
to 7800 = 27 mag, equivalent to 78i4 = 26.96 mag. They 
predict a total of 0.5 SN per WFPC2 field per visit. We 
therefore adopt an expected number of Np = 1 SN for every 
two-visit WFPC2 field reaching 78i4 = 27 mag. DF assume 
5.0 arcmin^ for the area of a WFPC2 field, while the value 
for the WF CCDs only is 4.69 arcmin^, so we scale their 
prediction by 0.938. 

In our calculations below, we will assume that 1/2 of 
the SNe with 7si4 < 27 mag have 78i4 < 26 mag, based 
on the results of S2000. We note that DF calculated the 
expected SN number counts in blank (non-cluster) fields, 
but because of the weak effect of cluster lensing on total 
background counts, these predictions can also be compared 
to our observations. 

Given a predicted SN surface density on the sky, we 
can calculate the expected number of SNe in the observed 
sample. For each data set, we determine the effective area S 
useful for SN search. This is just the sum of the areas that 
overlap the area of a reference epoch. For example, if useful 
data sets are available for two epochs with complete overlap, 
the effective area is just twice the illuminated area of the 
WF chips, or 9.38 arcmin^. In cases where one of the epochs 
cannot be searched, and was used only as a reference, its 
contribution to the effective area is zero. The predicted SN 
count per magnitude interval, dm, per observed cluster field 
is thus the surface density per magnitude interval, np{m)dm, 
times the fraction of a field effectively probed for SNe, times 
the search efficiency rj(m), with the latter modified by a 
colour term for data sets that were not observed in the 7 
band (see §3.3.2): 



n{m)dm = np(m) x 



9.38arcmin^ 



X rj{m) dm. 



(1) 



The total SN counts expected per field to a given limiting 

magnitude can be obtained by integrating Eq. 1 over m. 
S2000 and DF list their results for the surface density of 
SNe, Np{< rn), integrated up to some limiting magnitude, 
rather than the differential surface density np{m)dm. How- 
ever, since the SN counts are a steeply rising function of 
magnitude, we can obtain a good approximation of the to- 
tal by replacing the integral with the sum of two terms, for 
7814 < 26.4 (26) mag and for 26.4 < 78i4 < 27.4 (26 < 
7814 < 27) mag, according to the cumulative numbers 
given by S2000 (DF), and weighted by r?(26.4) [?7(26)] and 
77(27.4) [77(27)], respectively. The effective areas, efficiencies, 
and predicted SN numbers for each cluster field are given in 
Table 4. 

Summing over all fields, we obtain the total predicted 
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Table 4. Predicted Field SN Counts 



Cluster Eflfective area Efficiency Expected'^ SN number (A') 





(S) 




(ri) 






S2000 


S2000 


DP 


DP 






J — 26 


I = 26.4 


I = 27 


/ = 27.4 


/ < 26.4 


I < 27.4 


/ < 26 


/ < 27 


Abell 2218 


9.38 


0.48" 


0.47° 


0.46" 





0.58 


0.58 


0.23 


0.44 


Abell 1689 


7.81 


0.45 


0.29 








0.30 


0.30 


0.18 


0.18 


AC 114 


2.56 


0.69" 


0.62" 


0.5" 


0.41" 


0.21 


0.35 


0.09 


0.15 


MS1512.4+3647 


4.60 


0.62" 


0.53" 


0.37" 





0.32 


0.32 


0.14 


0.23 


MS2053.7-0449 


4.38 


0.77 


0.70 








0.40 


0.40 


0.17 


0.17 


MS1054.4-0321 


9.38 


0.84 


0.76 


0.65 





0.94 


0.94 


0.39 


0.70 


CL1604+4304 


7.94 


0.85 


0.75 


0.60 


0.56 


0.78 


1.37 


0.34 


0.58 


RXJ0848.6+4453 


9.38 


0.42^ 


0.31^ 








0.38 


0.38 


0.20 


0.20 


CL1645+46 


6.84 


0.82 


0.72 


0.51 





0.65 


0.65 


0.28 


0.45 



Notes: 

" Corrected with color term to convert to effective /-band efficiency (§3.3.2). 

* Corrected to account for loss of slowly-fading SNe in data with short time interval between epochs (§3.3.2). 
There are no offsets between sub-frames of this data set, making the verification of faint (/814 < 26.5) point 
sources impossible. The efficiency for (/814 > 26.5) is therefore set to zero. 

° The predictions of S2000 and DP, as given here, assume different cosmologies. Por inter-comparison, the 
predictions of S2000 should be scaled down by a factor of 0.7 (see text). 



SN numbers in the sample according to the models of S2000 
and DF. These predictions are compared with our observa- 
tions in Table 5. 

Looking at Table 5, one sees that the counts predicted 
by both DF and S2000 are formally consistent with the ob- 
servations. We note, however, that the predictions of S2000 
are always higher than the observed values. Moreover, the 
observed numbers given in Table 5 are upper limits, since, 
however unlikely, some of the apparent SNe we discuss may 
not be SNe (§ 3.2). Also, the "non-cluster" SN sample in- 
cludes objects (SNe 1996cp and 1996cq) for which the red- 
shift is uncertain, and hence may actually be cluster events. 
Finally, (§ 3.3.2) our adopted efficiency values are likely un- 
derestimates, leading to too-low predictions. For instance, if 
SN 1996cp is, in fact, a cluster event (§ 3.1.2), the differ- 
ence between the prediction of S2000 and the data becomes 
statistically significant. 

It is therefore instructive to investigate the sources of 
the differences between the predictions of S2000 and DF. 
The two groups assume very similar star-formation histories 
for the Universe (e.g., compare Figure 1 of DF with Figure 
2 of S2000). Each group assumes somewhat different forms 
for the delay functions of SNe-Ia (i.e., the SN-Ia rate as a 
function of time after an initial burst of star formation), 
but at these magnitudes, la's are a small fraction, 15-25 
per cent, of all SNe, so this has little effect on the total 
counts. For the UV spectrum of type-II-P SNe during their 
plateau phase, DF use a truncated blackbody model, while 
S2000 use an empirical spectrum, but their details are very 
similar, and in any case, this is of little consequence, since at 
/ < 27 mag the median redshift of the SNe is « ~ 0.7, where 
the I band probes rest-frame ~ 4800 A. The Q,m = 0.3, 
Q,A = 0.7 cosmology assumed by S2000 leads only to a ~ 30 
per cent increase in the total SN number to 7 < 27 mag, 
over that obtained using DF's assumed Q,m = 1 cosmology 
(P. Nugent, private communication). 

The bulk of the discrepancy between the predictions 
of S2000 and DF stems from two subtle differences in the 



treatment of type-II SNe, for which S2000 predict ~ 2 times 
more counts than DF. First, S2000 use a larger value for 
the dispersion of peak absolute magnitudes of type II-L SNe 
(cr = 1.3 mag), derived from the full sample of local type II- 
L SNe discussed by Miller & Branch (1990), and including 
the "overluminous" II-L SNe 1979C and 1980K. DF use a 
smaller dispersion (u = 0.5) derived from a sub-sample ex- 
cluding the two overluminous SNe, but retain the (brighter) 
peak luminosity derived from the full II-L SN sample. The 
large dispersion chosen by S2000 results in a high-luminosity 
tail of type II-L SNe, a large number of which will be visible 
in a flux-limited SN search. Second, S2000 adopted a con- 
stant dust extinction value for all type-II SNe, of Av = 0.45 
mag, equivalent to Ab = 0.66 mag. DF used a minimum 
extinction value oi Ab = 0.32 mag, lower than the one used 
by S2000, but modified it according to a random oriental 
tion of the host galaxy. This produced a large extinction 
for nearly edge-on hosts, with as many as 40 per cent of 
the core-collapse SNe suffering extinctions that are greater 
than 1 mag (T. Dahlen, private communication). This again, 
makes DF's SNe fainter than those of S2000, and the total 
predicted numbers lower. 

Returning to our observational results, it appears that 
the combined assumptions of DF for the type-II SN peak 
magnitudes, dispersions, and dust extinctions may be a bet- 
ter description of the high— SN population than those 
used by S2000. Of course, other variations on these assump- 
tions can be devised, which will produce the same predicted 
counts. Furthermore, some of the assumptions that are com- 
mon to the models of both DF and S2000 could also be 
wrong, but counterbalanced by wrong SN-II parameters. 
With improved knowledge of the relative fraction of differ- 
ent SN types, their luminosity functions, and their spectra, 
the observed SN counts will be able to constrain some of 
the other parameters, such as the extinction of core-collapse 
SNe at high redshift. 
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Table 5. Observed vs. Predicted Field SN Counts 



/814 Observed (total) Observed (non-cluster) Expected (S2000) Expected (DF) 
mag all bands /814 all bands /814 all bands /814 all bands /814 



< 26.4 


4 


3 


2 


2 


4.56 


3.45 






< 27.4 


5 


4 


3 


3 


5.28 


4.03 






< 26 


3 


2 


1 


1 






2.01 


1.55 


< 27 


5 


4 


3 


3 






3.09 


2.27 



3.4 Cluster supernova-la rates 

3.4- 1 Observed clusters and cluster supemovae 

In this section, we deal with the type-la SNe in the clusters 
themselves. As opposed to the preceding discussion of the 
field SN counts, which are dominated by type-IPs, the prop- 
erties of SNo-Ia are homogeneous and well characterised. As 
a result, wc can translate the observed numbers into ab- 
solute SN-Ia rates in rich clusters. We separate the clus- 
ters into two subsamples, a low-« sample consisting of four 
clusters at 0.18 < 2; < 0.37, and a high-z sample of three 
clusters at 0.83 < z < 1.27. Wc exclude from the calcula- 
tion the cluster MS2053.7— 0449, which is at an intermedi- 
ate redshift {z = 0.58), and whose data are relatively shal- 
low. The existence of the candidate cluster CL1645-I-46 at 
z > 1 was postulated by Jones et al. (1997), following the 
detection in this direction of a Sunyacv Zcldovich decre- 
ment, as well as a z = 3.8, 198-arcsec separation quasar 
pair, PC1643-I-4631A&B, which was suggested to be the 
lensed images of a single source. However, Kneissl, Sunyaev, 
& White (1998) argue that this cluster must be at 2; > 2.8, 
based on an improved upper limit on the X-ray emission. 
Since the observations in our sample cannot detect SNe at 
such high redshifts, this field is automatically excluded from 
our calculations of chister SN rates. 

The area of each cluster covered by the observations 
varies within our sample, with only the cluster core typi- 
cally covered for the low— z clusters, but also some of the 
cluster outskirts for the high— 2 clusters. We therefore limit 
our calculations of the SN rates to a metric radius of 500h^Q 
kpc around the cluster centres. This region is small enough 
that significant parts of it are covered by the WFPC2 im- 
ages available for the low— 2: clusters, but large enough to 
include much of the clusters' luminosity. SN 1996cp, which 
falls slightly outside this arbitrary region, is excluded from 
the calculations. 

Among the SNe we have found within the core regions, 
two SNe (1999gv and 1996cl) are most probably type-la's 
in their respective clusters, with one occurring in our low- 
z cluster sample and one in our high-2 sample. The high- 2; 
sample also includes SN 1996cq, which does not have a host 
with known redshift, and so may or may not be a cluster 
SN-Ia. If 1996cq belongs to the z — 0.83 cluster in whose 
field it was found, it is likely a la as well, since most core- 
collapse SNe would be as bright as this event only briefly at 
this redshift. We will therefore assume we have detected one 
cluster SN-Ia in the low-2 sample, and one or two cluster 
SNe-Ia in the liigli-z sample. The 95 per cent confidence 
ranges on the Poisson expectation values for the number of 
events in each sample, N, are then 0.07 < N < 4.75 at low 
z, and 0.07 < N < 6.3 at high z, where in the latter case we 



have taken the union of the 95 per cent probability ranges 
for observing i < 1 and i > 2 events. We propagate these 
uncertainties into our rate calculations, below. 

3.4-2 Rate calculation 

The SN rate per stellar luminosity unit in a sample of clus- 
ters is obtained by calculating 

As detailed below, N is the observed total number of cluster 
SNe detected. Ati is the effective visibility time of a SN- 
Ia at the cluster redshift, fi denotes the part of the rest- 
frame J3-band luminosity of a cluster, Lsi, that is probed 
by our SN search, and the summation is over all observed 
cluster epochs that are useful for a SN search. By convention, 
R is often expressed in SNu units, where 1 SNu = 1 SN 
century-i(lO^°Ls0)-\ 

The effective visibility time (sometimes called the "con- 
trol time") for a given epoch, Aii, is the period during which 
a SN-Ia is above the detection threshold of the data set, 
weighted by the detection efficiency as a function of magni- 
tude, r){m), during this period: 



(3) 



Here, m(t) is the SN light curve in the observed band, 
w.peafe,z is the maximum-light apparent magnitude of a SN 
at the redshift of the cluster, and mum is the limiting mag- 
nitude of the observation (?7 = for m > mum)- SN-Ia light 
curves in the specific WFPC2 filters and redshifts appropri- 
ate to our cluster sample were kindly provided by P. Nxigent, 
based on the data and calculations in Nugent et al. (2001). 
To allow numerical integration of Eq. 3, r]{m), which was 
estimated at discrete values of m in our detection simula- 
tions (§3.3.2), was linearly interpolated to a finer grid. The 
effective visibility times of the sample clusters are listed in 
Table 5. Note that, due to cosmological time dilation, a given 
time interval in the observer's frame corresponds to a time 
shorter by 1 -|- 2 in the rest frame of a cluster at z. However, 
each SN light curve is also stretched by 1 -|- 2;, and hence is 
detectable for a time that is 1 + 2 longer. Due to this can- 
cellation, 1-1-2 time dilation does not figure in A^, the total 
number of SNe detected. The rest-frame SN rate will there- 
fore be proportional to N divided by the rest-frame effective 
visibility time, as calculated by using, in Eq. 3, SN-la light 
curves without a cosmological 1 + z stretch in their time 
evolution. 

As is well documented (e.g., Hamuy et al. 1996; Riess 
et al. 1998; Perlmutter et al. 1997), SNe-Ia are not a per- 
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fcctly uniform population, but rather have a non-zero range 
of peak magnitudes and decline rates, with more luminous 
SNe fading more slowly and vice versa. Since we arc finding 
of order one or less SNe per cluster, this heterogeneity in the 
light curves of SNc-Ia introduces an uncertainty in the effec- 
tive visibility time. For example, in a given cluster, a faint, 
quickly declining, SN would have a shorter-than-average vis- 
ibility time, and one needs to account for the probability 
that such a SN is the one that was detected or missed in 
that cluster. To this end, we have created a family of SN-Ia 
light curves with "stretch" factors and correlated peak mag- 
nitudes. Based on the distribution observed by Perlmutter 
et al. (1999), wc assume SNc-Ia stretch factors between 0.8 
and 1.2, with a corresponding range in peak absolute mag- 
nitudes of —19.5 < Mb < —19.34. For each of the observed 
clusters, we then draw a light curve from the distribution 
and find the effective visibility time. The SN rate for the 
cluster sample is then calculated. This is repeated many 
times in a Monte-Carlo process, to gauge the effect of the 
light-curve non-uniformity on the derived SN rates. Wc find 
that the 95 per cent range of the effect on the final rates is 
about ±7 per cent for the low-z subsample, ±13 per cent for 
the higli-z subsample, and ±10 per cent for the entire sam- 
ple. Peculiar SNe-Ia (brighter, such as SN 1991T, or fainter, 
such as SN 1991bg; see, e.g, Li et al. 2001, and references 
within) are observed to be rare at high— a (Li et al. 2001), 
and so will have a negligible effect on our derived rates. 

The fraction, fi, of the luminosity probed at a given 
epoch is calculated by assuming that, within the 500h^Q kpc 
radius we consider, the surface brightness profile has a con- 
stant, 50h'^Q kpc core, and then falls as the projected radius 
to the —1 power. This is equivalent to assuming the light 
follows a cored isothermal mass distribution (e.g., Tyson & 
Fischer 1995). By integrating over the cluster area with use- 
ful coverage in the WFPC2 images, we obtain the fraction of 
the luminosity within 500/ij^,^ kpc that was actually searched 
for SNe. This fraction has a maximum value of / = 1 for 
a cluster in which both available epochs are useful for SN 
search, and the overlap area includes the entire 500/i5(/-kpc- 
radius region. Typical values of fi for our sample are 0.45 
for the low— z sub-sample, and 0.75 for the high— z clusters. 
The values of fi, multiplied by the useful number of epochs, 
for each cluster is given in Table 5. Reducing the adopted 
core radius by a factor of 5 changes the resulting SN rates 
only by about 5 per cent. Assuming a surface brightness pro- 
file that falls (unrcalistically) like the square of the radius 
lowers the SN rate by about 30 per cent. 

Normalizing the SN rate to the stellar luminosity re- 
quires knowledge of the cluster luminosities in the rest-frame 
B-band. We have found in the literature luminosities for 
several of the clusters in our sample. Where necessary, the 
published values have been corrected to our assumed cos- 
mology. We have taken the V-band luminosities for clus- 
ters with z ~ 0.3 and the /-band luminosities for clus- 
ters with z ~ 0.85 as reasonable representations of the rest 
frame -B-band luminosities. Our low— 2: sample consists of 
four rich clusters. Abell 2218 and AC114, which have pub- 
lished optical luminosities, are quite similar to Abell 1689 
and MS1512. 4+3647, which do not (this is also apparent 
from the //ST images used in this work), so we have assumed 
a value of Lb = 2.5 ±0.5 x 10^'^ h^^LBQ for the two clusters 
lacking published data, based on the Abell 2218 and AC114 



measurements. Our high— 2 sample includes two optically 
rich systems (MS1054.4-0321 and CL 1604-1-4304) shown to 
be at least as rich as the Coma cluster, with a mean value 
of Ls(< 500/i5o^kpc) ~ 5.5 x 10^^ Ko^bo- MS1054.4-0321 
does not have a published error estimate for its luminosity, 
and we assume an uncertainty of ±1.5 x 10^^ K^^Lbq - The 
third high-^; cluster, RXJ0848.6+4453 aX z = 1.27, presents 
a difficulty, as it is possibly poorer than the other two. How- 
ever, few examples of clusters at such high z exist, and esti- 
mates of typical luminosities have not been made. We there- 
fore adopt a plausible luminosity value for this cluster of half 
the mean luminosity of the other two, with a la uncertainty 
of one half the adopted value. This gives a probable range 
for the luminosity between zero (equivalent to dropping the 
cluster from the sample) and the high luminosity of the other 
two clusters. 

To assess the effect of the uncertainties of the cluster 
luminosities on the rate calculation, wo have included in the 
Monte-Carlo simulation described above, in addition to a 
light-curve draw, a luminosity draw for each cluster. Each 
cluster's luminosity was drawn from a Gaussian distribution 
with mean and a as given in Table 6. If the uncertainty in 
the visibility times is not included in the simulation (i.e., 
only a standard SN-Ia light curve is assumed) the effect of 
the luminosity uncertainty on the final rates is ±15 per cent 
at low z, +55 per cent, —30 per cent at high z, and +38 
per cent, —20 per cent for the entire sample (95 per cent 
confidence ranges). The large uncertainty in the luminosity 
of RXJ0848. 6+445 has a minor effect because, at its high 
redshift, the visibility time of SNe-Ia is short, giving it low 
weight in the calculation. The uncertainties in the luminosi- 
ties of the clusters clearly dominate over the uncertainties 
in the SN-Ia light-curve shape in the final SN-rate errors. 

With this input, we obtain from Eq. 2 a rest-frame clus- 
ter SN-Ia rate of /? = 0.20lo:?9ft50 SNu at 0.18 <z < 0.37 
and R = QAltl'^Mo SNu at 0.83 < z < 1.27, where the 
errors include the 95 per cent confidence Poissoii upper and 
lower limits, as well as the 68 per cent ranges due to the 
combined light-curve and luminosity uncertainties discussed 
above. Although the errors are large, we see no indication 
for any extreme evolution of the SN-Ia rate between the two 
redshift bins. We therefore also calculate a SN-Ia rate based 
on the two or three SNe in the combined cluster sample, with 
0.18 < z < 1.27, and a mean redshift (weighted by the rest- 
frame effective visibility time of each cluster) oi z = 0.41. 
We find R = 0.30lo:28/iio SNu, with the errors having the 
same meaning as above. 

3. 4.- 3 Comparison to models and to other measurements 

Table 7 summarises our rate measurements and compares 
them to those available in the literature. In local elliptical 
galaxies (which are the dominant population in the central 
regions of rich clusters), Cappellaro et al. (1997) find R = 
0. 067+0. 03/t5o SNu. Among distant field galaxies, measured 
rates are R = 0.20t";i9/iio SNu at z = 0.38 (Pain et al. 1997) 
and R = 0.16lg:^^/iio SNu at 2 = 0.55 (Pain et al. 2002; see 
Fabbro et al. 2000), where the quoted errors are la. The 
SN-Ia rates we have measured in the central regions of rich 
clusters at 0.18 < z < 1.27 are not markedly different from 
the rates in the local elliptical population, or in the general 
field population at similar redshifts. 
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Table 6. Cluster Luminosities and SN-Ia Visibility Times 



Cluster 




/iX 


At? 




[lOi^ft-fLBo] 


No. epochs'" 


[rest-frame days] 


Abell 2218 


2.5 ±0.5 


1.00 


219 


Abell 1689 


2.5 ±0.5 


0.80 


200 


AC 114 


2.5 ±0.25 " 


0.99 


208 


MS1512.4+3647 


2.5 ±0.5 


0.94 


147 


MS1054.4-0321 


5.5 ±1.5 f 


1.75 


101 


CL 1604+4304 


6.5 ±2 1 


1.00 


105 


RXJ0848.6+4453 


3 ±1.5 


1.50 


19 



Notes: 

" Rcst-framc B-band luminosity within R < 500h^Q kpc, and adopted la error, assuming f2m = 0.3, 
flA = 0.7. 

^ Fraction fi of the luminosity surveyed for SNe, times the number of epochs used for search. 

Effective visibility time, in days, of SNe-Ia at the cluster redshift, in cluster rest-frame. 

Squires et al. 1996 

Natarajan et al. 1998 
f Hoekstra, Pranx & Kuijken 2000 
9 Postman, Lubin, & Oke 2001 

No measurement available in the literature, see text. 



Table 7. SN-Ia Rates 



Environment 


Redshift 


SN Rate" 
/i§o[SNu] 


Reference 


Cluster (low— z) 
Cluster (high— z) 
Cluster (all) 


0.25 
0.90 
0.41 


20+"-^° 


This work 


E-SO 


local 


0.07 ±0.03 


CappcUaro 


SOa-Sb 


local 


0.09 ±0.03 


ct al. (1997) 


Sbc-Sd 


local 


0.11 ±0.04 




Cluster 


0.06 




Reiss (2000) 


Field 
Field 
Field 
Field 


0.11 

0.10 
0.38 
0.55 


10+"'"' 

11+009 

20+0 " 

16+°03 
^■'-"-0.02 


Hardin ct al. 2000 

Pain et al. 1997 
Fabbro et al. 2000 



Note: " All quoted errors are Icr. 

As mentioned in §1, a major motivation for measuring 
cluster SN rates is the possibility of addressing directly the 
puzzle of the large mass of iron and the excess entropy ob- 
served in clusters. For example, using gas-dynamical models 
for ICM enrichment by SNe, Brighenti & Mathews (1998) 
find that the cluster SN-Ia rate must be > 1.2/i?o SNu today 
and > 2.4/150 SNu at 2; = 1 to explain production of most 
of the iron in the ICM with SNe-Ia. These values are above 
our 95 per cent upper limits on the cluster SN-Ia rates. We 
note, however, that the uncertainty in our measured rates is 
dominated by the large Poisson errors resulting from small 
number statistics. Even a slightly larger cluster SN sample 
could more critically test the validity of the SN-Ia ICM en- 
richment scenario. 



4 CONCLUSIONS 

We have conducted a search for SNe using data from deep 
observations of galaxy clusters from the HST archive, and 



have detected six point-like transient events. Based on their 
observed characteristics and the known properties of differ- 
ent classes of variable objects, we have argued that most 
or all of these events are SNe. The redshifts and types of 
the likely host galaxies axe available in the literature for five 
events, demonstrating the benefits of observing well-studied 
fields. 

Folding in the observational parameters of our search, 
we have compared the numbers and magnitudes of the likely 
non-cluster (i.e., background or foreground) SNe we have 
found to published predictions of SN number counts. Such 
predictions depend strongly on a number of poorly known 
ingredients, and in particular the luminosity function and 
extinction of core-collapse SNe. The observed counts from 
our survey provide some initial constraints on these models. 
The observed counts, which are comparable to the predic- 
tions of DF, but somewhat lower than those of S2000, sug- 
gest that at z ~ 0.7, core-collapse SNe are relatively faint, 
due to cither extinction or to a luminosity distribution with- 
out a high-luminosity tail. However, being based on small 
SN numbers, these results should be treated with caution, 
and larger surveys, yielding significantly more events, are 
definitely called for. 

We have used the candidate type-la SNe in the clus- 
ters to estimate the SN-Ia rate in the central 500ftg(/ kpc 
of medium-redshift (0.18 < z < 0.37) and high-redshift 
(0.83 < z < 1.27) cluster sub-samples. The measured rates 
are consistent with SN-Ia measurements in field environ- 
ments, both locally and at high redshift. Our results ar- 
gue against an environmental dependence of the SN-Ia rate. 
Specifically, the observed low SN rate is in conflict with a 
scenario where SNe-Ia are the primary source of the large 
amounts of iron seen in the ICMs of massive clusters. 

This work shows that useful limits on high— z SN rates 
and counts can be derived from repeated, deep iJST imaging, 
even without any follow-up observations. The current dom- 
inant source of uncertainty is the Poisson error due to small 
number statistics. Next in importance are the lack of identi- 
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fication for some of the SNe, and the uncertainties in cluster 
luminosities. Obviously, similar analysis of a larger data set 
would yield SN rates and counts with higher precision, and 
planned observations (rather than an archival search) would 
enable some foUowup work that would reduce the systematic 
errors. Such an expanded program would also improve the 
prospects of discovering and utilising gravitationally lensed 
SNe, as outlined in §1. 

Excluding the search for SNe in the HDF by Gilliland 
et al. (1999), our search is the deepest SN survey con- 
ducted, and covers about 6 times the effective area of the 
HDF search. We find, on average, roughly one SN per deep 
(-^814 > 26 mag) duplicate WFPC2 cluster field, with about 
half of the SNe in the clusters and half in the field. At these 
magnitudes, there is a high density of SNe on the sky, de- 
tectable in significant numbers even in the small field of view 
of the HST cameras. Finding these SNe and characterising 
them, even if only in an incomplete and statistical way, can 
provide input and constraints to many current issues in as- 
tronomy and cosmology. If you are about to carry out deep 
imaging, whether from space or from the ground, do not 
do it all at once! Numerous SNe are out there, and can be 
found simply by splitting the observations into two or more 
well-separated epochs. 
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